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Abstract 


This  report  presents  a  library  of  wave-induced  ship  motion  modules  for  use  in  simulation 
applications.  The  library  is  intended  for  ship  defence  models  and  other  simulation  ap¬ 
plications  requiring  robustness,  fidelity,  and  high  computational  speed.  Ship  motions  in 
long-crested  random  seaways  arc  determined  using  ship  motion  response  amplitude  opera¬ 
tors  (RAOs)  computed  by  DRDC  Atlantic’s  SHIPM07  program.  The  output  ship  motions 
arc  given  as  perturbations  relative  to  the  steady  or  slowly  varying  input  ship  course.  Time 
domain  predictions  from  the  library  have  been  verified  through  comparisons  with  frequency 
domain  predictions  from  SHIPM07,  which  has  undergone  extensive  validation.  The  library 
is  available  as  a  Python  module  and  also  as  a  dynamically  linked  library  (DLL)  callable 
from  C++  and  other  compiled  languages.  During  benchmarking  on  a  800  MHz  Pentium 
III  computer,  ship  motion  predictions  required  approximately  0.002  s  per  time  step,  giving 
performance  50  times  faster  than  real  time  for  a  typical  simulation  application  using  a  time 
step  size  of  0.1  s. 


Resume 


Ce  rapport  presente  une  bibliotheque  de  modules  de  simulation  des  mouvements  des  na- 
vires  induits  par  des  vagues.  Cette  bibliotheque  a  ete  concuc  pour,  notamment,  des  ap¬ 
plications  de  simulation  de  la  defense  des  navires  qui  exigent  robustesse,  fidelite  et  une 
grande  rapidite  de  calcul.  Les  mouvements  des  navires  dans  une  mer  a  vagues  aleatoires 
de  grande  amplitude  sont  determines  au  moyen  d’operateurs  d'amplitude  de  reaction  des 
mouvements  du  navire  (RAO  pour  response  amplitude  operators)  calcules  au  moyen  du 
programme  SHIPM07  de  RDDC  Atlantique.  Les  mouvements  du  navire  produits  par  ce 
programme  sont  exprimes  sous  la  forme  d’ecarts  par  rapport  a  la  route  initiale,  qui  est 
constante  ou  varie  lentement.  Les  previsions  dans  le  domaine  temporel  obtenues  au  moyen 
de  la  bibliotheque  ont  ete  verifiees  en  effectuant  des  comparaisons  avec  des  previsions  dans 
le  domaine  des  frequences  etablies  avec  SHIPM07,  qui  a  ete  sounds  a  un  programme  de 
validation  complet.  Cette  bibliotheque  est  disponible  sous  la  forme  d’un  module  Python,  et 
egalement  sous  la  forme  d’une  DLL  (dynamically  linked  library  ou  bibliotheque  de  liens 
dynamiques),  qui  peut  etre  appelee  depuis  C++  et  d’autres  langages  compiles.  Lors  des 
essais  effectues  sur  un  ordinateur  Pentium  III  cadence  a  800  MHz,  les  previsions  des  mou¬ 
vements  du  navire  ont  necessite  environ  0,002  s  par  accroissement  de  temps,  ce  qui  est  50 
fois  plus  rapide  que  les  previsions  obtenues  en  temps  reel  avec  une  application  de  simula¬ 
tion  typique  fonctionnant  avec  un  accroissement  de  temps  de  0,1  s. 
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Executive  summary 

Introduction 

The  Canadian  Forces  Maritime  Warfare  Centre  tasked  DRDC  Atlantic  with  devel¬ 
oping  a  library  of  wave-induced  ship  motion  modules  for  use  in  simulation  appli¬ 
cations.  Fidelity,  robustness,  and  high  computational  speed  were  primary  consider¬ 
ations  when  selecting  a  suitable  approach. 


Principal  Results 

Ship  motions  in  long-crested  random  seaways  are  determined  using  ship  motion 
response  amplitude  operators  (RAOs)  computed  by  DRDC  Atlantic’s  SHIPM07 
program.  The  output  ship  motions  are  given  as  perturbations  relative  to  a  steady 
or  slowly  varying  input  ship  course.  Time  domain  predictions  from  the  library 
have  been  verified  through  comparisons  with  frequency  domain  predictions  from 
SHIPM07,  which  has  undergone  extensive  validation.  During  benchmarking  on 
a  800  MHz  Pentium  III  computer,  ship  motion  predictions  required  approximately 
0.002  s  per  time  step,  giving  performance  50  times  faster  than  real  time  for  a  typical 
simulation  application  using  a  time  step  size  of  0.1  s. 


Significance  of  Results 

The  library  will  bring  enhanced  fidelity  to  simulation  applications  for  which  ship 
motions  influence  system  performance.  Examples  of  entities  influenced  by  ship 
motions  include  sensors,  weapons,  ship-bome  aircraft,  and  humans.  The  fast  com¬ 
putational  speed  of  the  ship  motion  predictions  makes  it  practical  to  implement 
them  into  a  variety  of  simulation  applications. 


Future  Plans 

DRDC  Atlantic  is  developing  a  new  ship  motion  library  based  on  three-dimensional 
theory.  The  new  library  will  give  enhanced  fidelity  for  non-slender  hull  forms  and 
large  amplitude  motions.  The  new  library  will  also  be  able  to  simulate  ship  maneu¬ 
vering  in  response  to  input  rudder  and  propeller  states. 
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Introduction 

Le  Centre  de  guerre  navale  des  Forces  canadiennes  a  confie  a  RDDC  Atlantique  la  tache  de 
developper  une  bibliotheque  de  modules  de  simulation  des  mouvements  des  navires  induits 
par  des  vagues.  Dans  le  choix  de  la  methode,  la  fidelite,  la  robustesse  ainsi  qu’une  grande 
rapidite  de  calcul  constituaient  des  facteurs  determinants. 

Resultats  principaux 

Les  mouvements  des  navires  dans  une  mer  a  vagues  aleatoires  de  grande  amplitude  sont 
determines  au  moyen  d’operateurs  d'amplitude  de  reaction  des  mouvements  du  navire 
(RAO  pour  response  amplitude  operators)  calcules  au  moyen  du  programme  SHIPM07 
de  RDDC  Atlantique.  Les  mouvements  du  navire  produits  par-  ce  programme  sont  ex¬ 
primes  sous  la  forme  d’ecarts  par-  rapport  a  la  route  initiale,  qui  est  constante  ou  varie 
lentement.  Les  previsions  dans  le  domaine  temporel  obtenues  au  moyen  de  la  bibliotheque 
ont  ete  verifiees  en  effectuant  des  comparaisons  avec  des  previsions  dans  le  domaine  des 
frequences  etablies  avec  SHIPM07,  qui  a  ete  soumis  a  un  programme  de  validation  com- 
plet.  Lors  des  essais  effectues  sur  un  ordinateur  Pentium  III  cadence  a  800  MHz,  les  previsions 
des  mouvements  du  navire  ont  necessite  environ  0,002  s  par  accroissement  de  temps,  ce  qui 
est  50  fois  plus  rapide  que  les  previsions  obtenues  en  temps  reel  avec  une  application  de 
simulation  typique  fonctionnant  avec  un  accroissement  de  temps  de  0, 1  s. 

Importance  des  resultats 

Cette  bibliotheque  va  ameliorer  la  fidelite  des  applications  de  simulation  dans  lesquelles  les 
mouvements  des  navires  influent  sur  les  performances  du  systeme.  Citons,  a  titre  d’exemples 
d’entites  influencees  par  les  mouvements  des  navires,  les  capteurs,  les  armes,  les  aeronefs 
embarques  et  les  humains.  Grace  a  la  tres  grande  rapidite  d’etablissement  des  previsions 
concernant  les  mouvements  du  navire,  il  devient  possible  en  pratique  de  s’en  servir  dans 
toute  une  variete  d’ applications  de  simulation. 

Travaux  ulterieurs  prevus 

RDDC  Atlantique  est  en  train  de  developper  une  nouvelle  bibliotheque  de  simulations  des 
mouvements  des  navires,  basee  sur  la  theorie  tridimensionnelle.  Cette  nouvelle  bibliotheque 
va  ameliorer  la  fidelite  des  simulations  dans  le  cas  des  vagues  de  grande  amplitude. 
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1  Introduction 


Modelling  and  simulation  are  widely  used  for  a  variety  of  military  applications.  For 
systems  onboard  ships,  such  as  weapons  and  sensors,  wave-induced  ship  motions 
often  influence  system  effectiveness;  thus,  ship  motions  should  be  included  when 
modelling  ship  systems.  The  Canadian  Forces  Maritime  Warfare  Centre  tasked 
DRDC  Atlantic  with  developing  a  library  for  incorporating  ship  motions  into  ship 
defence  models. 

When  developing  a  suitable  approach  for  modelling  ship  motions,  the  main  consid¬ 
erations  were  robustness,  fidelity,  and  computational  speed.  The  selected  approach 
provides  ship  motions  in  long-crested  random  seaways  based  on  response  amplitude 
operators  (RAOs)  previously  computed  using  DRDC  Atlantic’s  SHIPM07  program 
[1].  SHIPM07  computes  ship  motions  in  the  frequency  domain  using  strip  theory, 
which  assumes  slender  hull  geometry.  Extensive  validation  of  SHIPM07  [2,  3, 4,  5] 
has  shown  very  good  agreement  with  model  tests  and  sea  trials  for  frigates  and  de¬ 
stroyers.  Good  results  can  generally  be  expected  from  SHIPM07  within  the  limi¬ 
tations  of  strip  theory.  Noticeable  deterioration  of  motion  predictions  can  occur  for 
non-slender  hull  forms  (e.g.,  length/beam  <  6)  or  at  higher  ship  speeds  (e.g.,  Froude 
numbers  >  0.4).  The  present  approach  assumes  linear  ship  motion  response,  which 
gives  good  results  for  frigates  in  seaways  up  to  NATO  sea  state  7. 

Section  2  of  this  report  presents  coordinate  systems  used  for  computations.  Sec¬ 
tion  3  describes  the  modelling  of  long-crested  random  seaways,  followed  by  Sec¬ 
tion  4  describing  the  modelling  of  ship  motions  in  those  seaways.  Section  5  presents 
the  software  implementation  of  the  model.  SHIPM07  predictions  of  response  am¬ 
plitude  operators  have  been  performed  for  three  ships  representative  of  a  corvette, 
frigate,  and  supply  ship,  as  described  in  Section  6.  Time  domain  predictions  have 
been  verified  with  frequency  domain  predictions  from  SHIPM07,  as  shown  in  Sec¬ 
tion  7.  The  main  report  finishes  with  conclusions  in  Section  8.  The  appendices  at 
the  end  of  the  report  are  intended  for  library  users  interested  in  importing  their  own 
ship  motion  RAOs. 

2  Coordinate  Systems 


The  coordinate  systems  for  the  present  library  were  selected  based  on  requirements 
from  the  Canadian  Forces  Maritime  Warfare  Centre.  Figure  1  shows  the  earth-fixed 
coordinate  system.  The  Xf  axes  is  positive  north  and  the  yj  axis  is  positive  east. 
Although  not  shown  in  Figure  1,  the  Zf  axis  has  its  origin  at  the  mean  still  water 
level  and  is  positive  down.  The  ship  heading  x  is  the  direction  to  which  the  ship 
is  going,  with  x  =  0  degrees  indicating  that  the  ship  is  heading  north.  Similarly, 
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Figure  1:  Fixed  Coordinate  System 

the  wave  direction  //  is  the  direction  to  which  the  waves  are  propagating,  with  ft  = 
0  degrees  indicating  that  the  waves  are  heading  north.  For  a  ship  travelling  with 
steady  speed  and  course,  the  nominal  ship  position  (not  including  perturbations  due 
to  wave-induced  motions)  will  be  as  follows: 

xf(t)  =  xf(t0)  +  (t  -  t0)  U  cos  x  (1) 

yf(t)  =  /(to)  +  (t  -  t0)  U  sin  x  (2) 

where  U  is  the  steady  ship  speed. 

Wave-induced  ship  motions  are  computed  as  perturbations  relative  to  a  steady  or 
slowly  varying  ship  course.  The  wave-induced  ship  motions  are  computed  using 
a  translating  earth  coordinate  system  (Figure  2)  that  moves  with  the  nominal  ship 
course.  Wave-induced  ship  motions  are  computed  for  all  six  degrees  of  freedom. 

As  mentioned  earlier,  wave-induced  ship  motions  are  evaluated  using  RAOs  pre¬ 
viously  computed  by  the  frequency  domain  program  SHIPM07.  RAOs  for  each 
motion  mode  are  functions  of  ship  speed  U,  relative  sea  direction  (3,  and  incident 
wave  frequency  uj.  Figure  3  shows  the  definition  of  relative  sea  direction,  which  is 
related  to  ship  heading  and  wave  heading  as  follows: 

P  =  F  ~  X  (3) 
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Figure  2:  Translating  Earth  Coordinate  System 


Due  to  symmetry,  ship  motion  RAOs  only  need  to  be  stored  for  relative  sea  direc¬ 
tions  of  0  to  180  degrees. 

3  Modelling  of  Long  Crested  Random 

Seaways _ 

The  modelling  of  seaways  influencing  ship  motions  is  discussed  in  detail  in  Ref¬ 
erence  6.  The  present  work  uses  the  coordinate  system  presented  in  the  previous 
section. 

The  surface  elevation  of  a  regular  seaway  can  be  modelled  as  follows: 

(i(xf  ,yf  ,t)  =  a  cos  kj  {x*  cos  v  +  i/  sin  //)  —  ujj  t,  —  ef  (4) 

where  a  is  the  wave  amplitude,  ki  is  the  incident  wavenumber,  x j  and  r//  are  the 
horizontal  plane  coordinates,  toj  is  the  incident  wave  frequency,  and  is  the  phase 
lead  angle  for  the  wave  trough  at  the  origin.  The  surface  elevation  uses  the  same 
sign  convention  as  the  z*  axis,  with  downward  being  positive.  The  present  work 
assumes  deep  water  (i.e.,  water  depth  greater  than  half  of  incident  wavelength),  and 
the  following  dispersion  relation  applies: 

ki  —  —  (5) 
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Figure  3:  Sea  Direction  Relative  to  Forward  Ship  Speed 
where  g  is  gravitational  acceleration. 

A  random  seaway  can  be  simulated  by  linear  superposition  of  a  finite  number  of 
regular  wave  components.  Chakrabarti  [7]  provides  a  useful  overview  of  this  ap¬ 
proach.  Using  the  fixed  axis  systems  of  Figure  1,  the  Zf  value  of  the  free  surface 
for  a  unidirectional  random  seaway  can  be  simulated  as  follows: 

Nj 

C i(xf,yf,t )  =  cos 

i= 1 

where  Nj  is  the  number  of  regular  wave  components  used  to  represent  the  random 
seaway,  a*  is  the  wave  amplitude  for  component  i,  ki-i  is  the  wavenumber  for 
component  i,  c is  the  incident  wave  frequency  for  component  i,  and  e{_i  is  the 
phase  for  wave  component  i.  The  wave  amplitude  for  component  i  is  based  on  a 
specified  wave  spectrum  as  follows: 

a*  =  \/2  (cc/_;)  A(c o^i)  (7) 

where  SWI  is  the  spectral  density  and  A(cc/_j)  is  the  wave  frequency  interval  for 
component  i.  The  wave  phase  components  e{_{  are  obtained  using  random  num¬ 
ber  generation  based  on  a  uniform  distribution  between  0  and  2tt.  The  following 
relationship  exists  between  the  the  wave  component  amplitudes  and  the  standard 
deviation  of  water  elevation: 

1  N< 

ff2(C)  =  2  £a?2  (8) 

i— 1 

The  fidelity  of  the  simulated  seaway  will  increase  with  number  of  wave  compo¬ 
nents.  A  minimum  of  19  components  is  typically  used  for  simulating  unidirectional 
random  seaways. 


(x*  cos  v  +  y*  sin  v)  —  t 


-I-i 


(6) 
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Reference  6  describes  several  commonly  used  wave  spectral  models.  The  Bretschnei- 
der  spectrum  is  the  most  widely  used  spectrum  for  waves  in  the  open  ocean,  and  has 
been  selected  for  the  present  application.  Based  on  the  15th  International  Towing 
Tank  Conference  (ITTC)  [8],  the  formulation  for  the  Bretschneider  spectrum  is: 


sUI  M 


486.0  H2S 

T4  cc5 

p 


exp 


—  1948.2’ 

.w. 


(9) 


The  above  spectrum  is  defined  in  terms  of  peak  wave  period  Tp,  which  is  associated 
with  the  wave  frequency  at  which  SUl(u>i)  has  its  maximum.  For  a  Bretschnei¬ 
der  spectrum,  the  following  relations  exist  between  the  peak  wave  period  and  the 
average  and  zero-crossing  wave  periods: 


T,  =  0.773  Tp  (10) 

Tz  =  0.710  Tp  (11) 


The  current  model  of  ship  motions  in  waves  requires  the  user  to  select  suitable 
values  of  significant  wave  height  Hs  and  peak  wave  period  Tp.  Available  wave 
climate  data  can  provide  useful  guidance  for  selecting  seaway  parameters.  Table  1 
gives  NATO  sea  state  numbers  from  Lee  and  Bales  [9]  and  associated  wave  and 
wind  parameters. 

Table  1:  Sea  State  Parameters  for  the  North  Atlantic  (reproduced  from 
Reference  9) 


Sea 

State 

Number 

Significant  Wave 
Height  (m) 

Sustained  Wind 

Speed  (knots) 

Percentage 

Probability 

of  Seastate 

Peak  Wave  Period 

Range 

Most 

Probable 

Range 

Mean 

Range 

Mean 

0-1 

0-0.1 

0.05 

0-6 

3 

0.7 

- 

- 

2 

0.1  -0.5 

0.3 

7  -  10 

8.5 

6.8 

3.3-12.8 

7.5 

3 

0.5  -  1.25 

0.88 

11  -  16 

13.5 

23.7 

5.0-14.8 

7.5 

4 

1.25-2.5 

1.88 

17-21 

19 

27.8 

6.1  -  15.2 

8.8 

5 

2.5-4 

3.25 

22-27 

24.5 

20.6 

8.3-15.5 

9.7 

6 

4-6 

5 

28-47 

37.5 

13.1 

9.8-16.2 

12.4 

7 

6-9 

7.5 

48-55 

51.5 

6.1 

11.8-  18.5 

15.0 

8 

9-  14 

11.5 

56-63 

59.5 

1.1 

14.2-  18.6 

16.4 

>8 

>  14 

>  14 

>63 

>63 

0.05 

18.0-23.7 

20.0 

When  specifying  input  values  for  significant  wave  height  and  peak  wave  period,  it  is 
important  that  the  resulting  seaway  have  realistic  wave  steepnesses.  Analysis  of  ob¬ 
served  wave  data  by  Buckley  [10]  indicates  that  the  following  condition  should  be 
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met  when  selecting  input  values  for  significant  wave  height  and  peak  wave  period: 


Hs 

Tl  g/( 2  tt) 


<  0.049 


(12) 


The  left  hand  side  of  the  above  equation  represents  a  nominal  wave  steepness,  with 
the  denominator  representing  the  wavelength  in  deep  water  associated  with  the  peak 
wave  period  Tp.  The  existence  of  steeper  waves  is  restricted  by  the  occurrence  of 
wave  breaking. 


4  Modelling  of  Ship  Motions  in  Long 
Crested  Random  Seaways 


When  computing  ship  motions  for  use  in  ship  defence  models,  it  is  important  that 
computations  can  proceed  much  faster  than  real  time  for  the  following  reasons: 

•  ship  defence  models  consist  of  many  different  elements  which  must  share  com¬ 
putational  resources, 

•  multiple  runs  of  scenarios  are  often  required  for  Monte  Carlo  simulation, 

•  many  different  scenarios  are  often  examined  to  determine  optimal  tactics. 

To  achieve  a  suitable  balance  of  fidelity  and  computational  speed,  ship  motions  in 
the  time  domain  can  be  evaluated  based  on  motion  response  amplitude  operators 
(RAOs)  previously  computed  in  the  frequency  domain.  For  the  present  applica¬ 
tion,  DRDC  Atlantic’s  SHIPM07  program  [1]  can  be  used  to  determine  frequency 
domain  RAOs.  Comparisons  with  model  tests  [4]  and  full  scale  trials  [5]  indicate 
that  SHIPM07  gives  very  good  predictions  for  naval  frigates  operating  in  seas  with 
significant  wave  heights  up  to  7  m. 

SHIPM07  produces  complex  response  amplitude  operators  of  wave-induced  mo¬ 
tions  for  a  ship  travelling  with  steady  speed  and  heading.  For  a  ship  in  a  sinusoidal 
seaway,  the  relationship  between  ship  motions  in  the  time  and  frequency  domains 
can  be  expressed  by: 

rjj(t)  =  Real  {a  r)*(U,  f3,u>i)  exp(iu;ef)}  for  j  =  1  —  6  (13) 

where  r/j  is  the  displacement  for  mode  j,  Real{. . .}  denotes  the  real  part  of  the 
complex  term  in  brackets,  and  ?/*  is  the  complex  motion  RAO  (i.e.,  the  complex 
motion  amplitude  in  unit  amplitude  waves)  for  mode  j.  Assuming  linear  motion 
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response,  the  complex  motion  RAO  is  a  function  of  ship  speed  U,  relative  sea  di¬ 
rection  (5,  and  incident  wave  frequency  u>i.  The  wave  encounter  frequency  ue  for 
the  ship  travelling  at  steady  speed  and  heading  is  given  by: 

cue  =  \u>i  —  U  kj  cos  /3\  (14) 


Although  Equation  (13)  assumes  linear  motion  response,  roll  response  is  actually 
nonlinear  due  to  the  importance  of  nonlinear  viscous  roll  damping.  To  account  for 
this  effect,  ship  motion  RAOs  can  be  computed  for  a  specified  nominal  seaway.  For 
frequency  domain  computations,  an  iterative  procedure  is  typically  used  to  obtain 
the  correct  roll  amplitude  in  a  specified  seaway. 


When  solving  for  ship  motions  in  the  frequency  domain,  computations  are  per¬ 
formed  in  translating  earth  axes  which  travel  at  the  steady  speed  and  heading  of  the 
ship.  In  the  present  coordinate  system,  it  is  specified  that  the  wave  trough  will  be 
located  at  x  —  0,  y  =  0  (see  Figure  2)  at  time  t  —  0.  The  resulting  motions  in  the 
time  domain  will  be  related  to  the  predicted  frequency  domain  motions  as  follows: 

rjj(t)  =  a  \Vj(U,  P,u>1)\  cos  [uje  t  +ej(U,(3,u;I)\  forj  =  l-6  (15) 


The  phase  lead  of  the  motions  relative  to  the  wave  trough  being  at  the  origin  is 
given  by: 


Imag  {r)j(U,  fi,u>i)} 

arctan - ^ - -r- 

Real  {ri*(UJ3,(jjr)} 


(16) 


where  Imag  {. . .}  denotes  the  imaginary  part  of  the  complex  term  in  brackets. 


When  computing  ship  motions  within  ship  defence  models,  it  is  essential  to  account 
for  the  possibility  of  ship  speed  and/or  heading  changing  with  time,  as  can  occur 
while  a  ship  maneuvers.  For  a  ship  in  a  fixed  seaway  modelled  using  Equation  (4), 
the  instantaneous  ship  motions  can  be  expressed  as: 

Nj 

Vj(t )  =  hi (£4/3,^/) I  COS  [u;e_i(f7,/?,u;/_j)  t  +  £/_,  +  ej\  (17) 

i— 1 

where  e/_j  is  the  phase  lead  angle  of  the  incident  wave  trough  of  wave  component 
i  in  translating  earth  axes  with  their  origin  at  the  ship  centre  of  gravity.  The  wave 
elevation  in  translating  earth  axes  from  wave  component  i  is  given  by: 

Cj -i(x,y,t)  =  a.i  cos  (x cos  f3  +  y sin/3)  -  cce_;  t  -  e/_j]  (18) 

Based  on  Equations  (6)  and  (18),  the  wave  elevation  phase  in  translating  earth  co¬ 
ordinates  at  any  time  during  a  simulation  is  given  by: 

ei-i  =  efj_t  +  t  (ui-i  -  uJe-i)  -  kT-i  (xf  cos  /i  +  yf  sin/x)  (19) 
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When  applying  the  above  equation,  xJ  are  yf  are  the  instantaneous  coordinates  of 
the  ship  centre  of  gravity.  If  both  the  ship  speed  and  heading  remain  constant  during 
a  simulation,  then  the  wave  elevation  phase  term  f /_,  will  also  remain  constant. 

Wave-induced  ship  velocities  and  accelerations  are  often  required  during  simula¬ 
tions.  These  terms  can  be  obtained  from  differentiation  of  Equation  (17): 

N ! 

7 j(t)  =  -Ve-i  ai  \v*j(U,(3,ui)\  sin  [ue-i(U,  (3,  +  e/_j  +  ej]  (20) 

i= 1 
Nj 

Vj(t)  =  |r7* (E/, /3, cu7) |  cos[cje_i(?7,/3,a;i_i)t  +  e/-»  +  e,-]  (21) 

1=1 

The  equations  presented  thus  far  give  the  motions  at  the  ship  centre  of  gravity.  Local 
motions  at  a  ship  location  offset  from  the  centre  of  gravity  are  often  of  interest.  For 
a  location  on  a  ship  specified  using  the  translating  earth  coordinates  of  Figure  2,  the 
local  motions  in  translating  earth  coordinates  are  given  as  follows: 


6(f) 

=  7i(f)  +  z  '75(f)  -  7  7e(f) 

(22) 

6(f) 

=  rj2 (t)  -  274(f)  +  xrj6(t) 

(23) 

6(f) 

=  73(f)  +  7  74(f)  -  x  75(f) 

(24) 

The  above  equations  are  based  on  the  assumption  of  small  angular  motions. 

5  Software  Implementation 


A  module  has  been  developed  using  the  computer  language  Python  [1 1,  12]  to  pro¬ 
vide  ship  motions  for  ship  defence  simulations.  Python  is  an  open  source  program¬ 
ming  language  that  has  a  wide  user  base,  including  many  scientific  users.  Python 
is  object-oriented  like  C++  and  Java;  however,  it  is  a  higher  level  language,  which 
results  in  significantly  fewer  lines  of  code  being  required  to  develop  a  specific  ap¬ 
plication.  The  availability  of  the  Numeric  library  [13]  has  been  a  primary  reason 
for  Python’s  acceptance  among  scientific  users.  Python  has  been  used  since  2001 
for  development  of  ship  motion  predictions  at  DRDC  Atlantic  [6,  14,  15]. 

For  the  present  application,  ship  motion  RAOs  are  computed  using  DRDC  At¬ 
lantic’s  SHIPM07  [1],  SHIPM07  is  based  on  strip  theory,  which  models  the 
geometry  of  a  hull  using  a  number  of  two-dimensional  sections  (typically  21)  along 
the  hull  length.  Strip  theory  is  based  on  the  assumptions  of  slender  hull  geom¬ 
etry  and  moderate  ship  speed.  Comparisons  with  experiments  and  full-scale  tri¬ 
als  [2,  4,  5,  16]  indicate  that  SHIPM07  gives  very  good  results  for  ships  with 
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length/beam  ratios  greater  than  6  travelling  at  Froude  numbers  less  than  0.4  (27 
knots  for  a  frigate  with  waterline  length  of  120  m). 

A  Python  class  ShipMotionDefenceRAODBSeaway  has  been  developed  to  provide 
ship  motions  for  ship  defence  models.  An  instance  of  this  class  is  initialized  with 
the  following  input  variables: 

•  fileNameSm7RAO  —  file  name  of  SHIPM07  ASCII  post-processing  file  with 
motion  RAOs, 

•  hSig  —  significant  wave  height  (m), 

•  tPeak  —  peak  wave  period  (s), 

•  waveHeadingToDeg  —  direction  (deg)  to  which  waves  are  heading  (0  degrees 
for  waves  heading  north,  90  degrees  for  waves  heading  east), 

•  waveFreqs  —  sequence  of  wave  frequencies  (rad/s).  If  these  are  not  given  as 
input,  then  default  values  of  0.2,  0.3,  . . .,  2.0  rad/s  are  used. 

•  wavePhaseSeeds  —  sequence  of  two  integer  seed  numbers  for  generation  of 
phases  for  seaway  components.  If  these  are  not  given  as  input,  then  default 
values  of  1001  and  7001  are  used. 

In  the  above  variable  descriptions,  the  term  “sequence”  indicates  a  collection  of 
values  which  can  be  in  the  form  of  a  Python  list,  tuple,  or  Numeric  array. 

Annex  A  gives  the  format  of  the  SHIPM07  post-processing  file.  When  preparing 
a  SHIPM07  post-processing  file  to  be  used  for  ship  motion  simulations,  attention 
must  be  given  to  selecting  suitable  ranges  of  ship  speeds,  relative  sea  directions,  and 
wave  frequencies.  The  SHIPM07  computations  must  encompass  the  range  of  ship 
speeds  that  will  be  required  for  motion  simulations.  It  is  recommended  that  a  ship 
speed  interval  of  5  knots  or  smaller  be  used.  SHIPM07  relative  sea  directions  must 
have  a  minimum  value  of  0  degrees,  a  maximum  value  of  180  degrees,  and  a  direc¬ 
tion  interval  no  greater  than  30  degrees.  It  is  recommend  that  wave  frequencies  of 
0.2,  0.3, . . .,  2.0  rad/s  be  used.  If  a  simulation  requires  motions  at  lower  frequencies 
than  are  available  in  the  SHIPM07  file,  then  approximations  are  made  based  on  the 
lowest  frequency  for  which  SHIPM07  computations  are  available.  If  a  simulation 
requires  motions  at  higher  frequencies  than  are  available  in  the  SHIPM07  file,  then 
ship  motions  are  assumed  to  be  zero. 

During  initialization  of  a  ShipMotionDefenceRAODBSeaway  object,  the  ship  mo¬ 
tion  RAOs  are  converted  from  the  SHIPM07  axis  system  of  Reference  1  to  that 
given  in  Figure  2.  The  conversion  process  must  account  for  the  convention  in  the 
present  work  that  motion  phases  are  relative  to  when  a  wave  trough  passes  the 
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Table  2:  Multiplication  Factors  for  Obtaining  Complex  Ship  Motion  RAOs  from 
SHIPM07  RAOs 


Surge 

-1 

Sway 

+1 

Heave 

+1 

Roll 

-1 

Pitch 

+1 

Yaw 

+1 

ship  CG,  with  conversion  factors  given  in  Table  2.  Motions  RAOs  are  ultimately 
stored  at  the  wave  frequencies  for  the  simulated  seaway.  Linear  interpolation  is 
used  when  wave  frequencies  from  the  SHIPM07  RAOs  differ  from  those  for  the 
simulated  seaway. 

Once  the  ShipMotionDefenceRAODBSeawoy  object  has  been  initialized,  motion 
displacements  at  the  ship  centre  of  gravity  at  a  given  instant  in  time  can  be  obtained 
using  the  method  displacementsCG,  which  has  the  following  input  variables: 

•  t  —  time  (s), 

•  speed  —  ship  speed  at  time  t  (m/s), 

•  shipHeadingToDeg  —  ship  heading  (deg)  (0  degrees  for  ship  heading  north,  90 
degrees  for  ship  heading  east), 

•  xf  —  x  coordinate  of  ship  CG  in  fixed  coordinate  system, 

•  yf  —  V  coordinate  of  ship  CG  in  fixed  coordinate  system, 

Method  displacementsCG  returns  the  following  wave-induced  motion  components 
in  the  translating  axis  system  of  Figure  2: 

•  surge  (m), 

•  sway  (m), 

•  heave  (m), 

•  roll  (deg), 

•  pitch  (deg), 

•  yaw  (deg). 
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Note  that  the  above  motion  displacements  are  relative  to  the  intended  ship  course 
in  calm  water.  Methods  velocitiesCG  and  accelerationsCG  are  also  available  for 
obtaining  wave-induced  ship  velocities  and  accelerations.  Velocities  have  units  of 
m/s  and  deg/s,  and  accelerations  have  units  of  m/s2  and  deg/s2. 

Benchmark  tests  have  been  performed  on  a  800  MHz  Pentium  III  computer  to  deter¬ 
mine  computational  requirements  for  modelling  ship  motions.  Each  call  to  method 
displacementsCG  requires  approximately  0.002  s  of  CPU  time.  For  a  representa¬ 
tive  simulation  computing  ship  motions  at  a  time  interval  of  0. 1  s,  the  ship  motion 
predictions  would  run  approximately  50  times  faster  than  real  time. 

To  facilitate  usage  by  programs  written  in  C++,  the  Python  class  ShipMotionDe- 
fenceRAODBSeaway  has  been  implemented  as  a  C++  class  in  a  dynamically  linked 
library  (DLL).  The  DLL  interface  is  given  in  the  C++  header  file  DRDCShipMo- 
tionDefence.h,  and  an  example  C++  program  testShipMotionDefence.exe  demon¬ 
strates  usage.  Execution  of  the  ShipMotionDefenceRAODBSeaway  DLL  and  test 
program  requires  that  the  files  boost_python.dll,  DRDCPythonWrapper.dll,  and  stl- 
port_vc646.dll  be  available  in  the  current  path  at  runtime. 

6  Sample  Ships  for  Pre-Computed  Ship 
Motion  RAOs 


To  facilitate  integration  of  ship  motions  into  ship  defence  models,  files  with  motion 
RAOs  have  been  prepared  using  SHIPM07  for  three  sample  ships  of  different  sizes 
with  properties  given  in  Table  3.  The  three  ships  are  geosims,  with  the  120  m  frigate 
being  the  parent  ship  geometry,  and  the  70  m  corvette  and  the  200  m  supply  ship 
being  scaled  from  the  frigate  dimensions.  The  hull  geometry  is  taken  from  Refer¬ 
ence  17,  and  is  representative  of  modem  frigates,  which  have  wide  transom  stems. 
Figure  4  shows  the  hull  lines  and  waterline.  Each  sample  ship  has  bilge  keels,  a 
skeg,  and  a  single  rudder.  Appendage  dimensions  were  chosen  to  be  representative 
of  modem  naval  vessels.  The  height  of  the  center  of  gravity  for  each  ship  is  also 
representative  of  modem  naval  vessels,  and  is  primarily  influenced  by  requirements 
of  stability  standards  for  prevention  of  ship  capsize.  In  summary,  the  hull  forms, 
appendages,  and  inertial  properties  of  the  sample  ships  are  representative  of  op¬ 
erational  naval  vessels,  and  the  resulting  motions  in  waves  can  be  expected  to  be 
realistic.  Annex  B  gives  the  SHIPM07  input  file  for  the  sample  frigate. 

The  computed  ship  motion  RAOs  for  the  sample  ships  cover  the  following  parame¬ 
ter  ranges: 

•  ship  speeds  U  of  0,  5,  10,  . . .,  30  knots, 
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Table  3:  Properties 

of  Sample 

Ships 

Corvette 

Frigate 

Supply  Ship 

Length  between  perpendiculars  L  (m) 

70.0 

120.0 

200.0 

Beam  B  (m) 

8.228 

14.106 

23.513 

Draft  at  midships  Tmid  (m) 

2.450 

4.200 

7.0 

Trim  by  stem  ts  (m) 

0.000 

0.000 

0.000 

Displacement  A  (tonnes) 

741 

3735 

17295 

Height  of  CG  above  waterline  (m) 

1.050 

1.800 

3.0 

Metacentric  height  GM  (m) 

0.815 

1.398 

2.334 

Roll  radius  of  gyration  (m) 

2.800 

4.800 

8.000 

Pitch  radius  of  gyration  (m) 

17.5 

30.0 

50.0 

Yaw  radius  of  gyration  (m) 

17.5 

30.0 

50.0 

Natural  roll  period  (s) 

7.1 

9.1 

11.7 
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•  relative  sea  directions  (3  of  0,  15,  30,  . . 180  degrees, 

•  incident  wave  frequencies  c uj  of  0.2,  0.3,  0.4,  . . .,  2.0  rad/s. 

The  ship  speeds  are  considered  to  cover  a  realistic  range  of  values.  The  relative 
sea  directions  were  selected  such  that  all  relative  headings  are  covered  due  to  the 
assumption  of  lateral  ship  symmetry.  The  incident  wave  frequencies  encompass  all 
relevant  frequencies  at  which  significant  energy  can  occur  to  influence  ship  mo¬ 
tions. 


7  Verification  with  Frequency  Domain 
Motion  Computations 


As  indicated  previously,  the  ship  motion  program  SHIPM07  has  been  validated 
extensively,  including  comparisons  with  model  tests  [2,  4]  and  full-scale  trials  [5]. 
These  comparisons  indicate  that  SHIPM07  gives  good  agreement  for  slender  hull 
forms  travelling  at  moderate  forward  speed. 

To  verify  the  present  library  for  ship  motions  in  the  time  domain,  comparisons 
of  motions  have  been  made  with  frequency  domain  predictions  from  SHIPM07. 
Figures  5  and  6  show  motions  at  the  ship  centre  of  gravity  in  NATO  sea  state  5,  and 
indicate  excellent  agreement.  To  verify  that  the  relative  phasing  of  motion  modes 
is  correct,  further  computations  have  been  performed  for  a  point  offset  from  the 
ship  centre  of  gravity.  The  selected  point  is  located  at  station  3  (stations  0  and  20 
represent  the  forward  and  aft  perpendiculars  respectively),  has  a  lateral  offset  of  5  m 
to  port,  and  is  2  m  above  the  ship  baseline.  The  lateral  and  vertical  displacements  in 
Figures  7  and  8  show  excellent  agreement  between  the  frequency  and  time  domain 
predictions,  indicating  that  phasing  between  motion  modes  is  correct.  To  verify 
phasing  between  the  simulated  ship  motions  and  seaway,  relative  motions  (local 
vertical  displacement  minus  local  seaway  elevation)  have  been  computed  for  the 
offset  location  used  for  verifying  lateral  displacement.  Figure  9  shows  excellent 
agreement  between  the  results  from  the  frequency  and  time  domains. 
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RMS  heave  (m)  RMS  sway  (m) 


Time  domain 
Frequency  domain 


Figure  5:  RMS  Motions  and  Zero  Crossing  Periods  for  Frigate  Surge,  Sway,  and 
Heave,  Hs  =  3.25  m,  Tp  =  9.7 
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RMS  pitch  (deg) 


-  Time  domain 


Frequency  domain 


Relative  sea  direction  (3  (deg)  Relative  sea  direction  /3  (deg) 


Relative  sea  direction  [3  (deg)  Relative  sea  direction  (3  (deg) 


Relative  sea  direction  (3  (deg)  Relative  sea  direction  (3  (deg) 


Figure  6:  RMS  Motions  and  Zero  Crossing  Periods  for  Frigate  Roll,  Pitch,  and  Yaw, 
Hs  =  3.25  m,  Tp  =  9.7 
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Time  domain 
Frequency  domain 


Figure  7:  RMS  Motions  and  Zero  Crossing  Periods  for  Lateral  Displacement  at 
Offset  Point  on  Frigate,  Hs  =  3.25  m,  Tp  =  9.7 


Time  domain 
Frequency  domain 


Figure  8:  RMS  Motions  and  Zero  Crossing  Periods  for  Vertical  Displacement  at 
Offset  Point  on  Frigate,  Hs  =  3.25  m,  Tp  =  9.7 
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Time  domain 
Frequency  domain 


Figure  9:  RMS  Motions  and  Zero  Crossing  Periods  for  Relative  Vertical 
Displacement  at  Offset  Point  on  Frigate,  Hs  =  3.25  m,  Tp  =  9.7 
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8  Conclusions 


A  library  has  been  developed  for  predicting  wave-induced  ship  motions  in  long- 
crested  seas.  The  library  is  suitable  for  ship  defence  models  and  other  applications 
requiring  a  combination  of  good  fidelity  and  high  computational  speed.  Wave- 
induced  motions  are  computed  as  perturbations  from  a  prescribed  ship  trajectory. 

Using  a  relatively  slow  desktop  computer,  computations  of  ship  motions  require  ap¬ 
proximately  0.002  s  of  CPU  time  per  time  step.  For  a  typical  simulation  application 
with  a  step  size  of  0.1  s,  the  ship  motion  computations  will  run  approximately  50 
times  faster  than  real  time. 

To  facilitate  usage  of  the  library,  ship  motion  RAOs  are  available  for  three  repre¬ 
sentative  ships  with  lengths  of  70  m  (nominal  corvette),  120  m  (nominal  frigate), 
and  200  m  (nominal  supply  ship). 
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Symbols  and  Abbreviations 


a 

Qi 

B 

CG 

DLL 

Fn 

GM 

9 

Hs 

ki 

ki-i 

L 

RMS 

sUI  M 

7" 1 

-*■  mid 

T 

1p 

Tz 

T\ 


to 

U 

x,y,z 

xf,yf,zf 

(3 

e3 

£/ 

ef 

ti 

J 

ti~~i 

C / 

Vj 

$ 

9 

Zj 

a 

X 

Ui 

Wl-i 


wave  amplitude 

amplitude  of  wave  component  i 
beam 

centre  of  gravity 

dynamically  linked  library 

Froude  number 

metacentric  height 

gravitational  acceleration 

significant  wave  height 

incident  wavenumber 

incident  wavenumber  of  wave  component  i 

ship  length  between  perpendiculars 

root-mean- square 

wave  spectral  density 

draft  at  midships 

peak  wave  period 

zero-crossing  period  of  spectrum 

average  wave  period 

trim  by  stern 

reference  time 

ship  forward  speed 

coordinates  in  translating  earth  axes 

coordinates  in  earth-fixed  axes 

wave  direction  relative  to  ship 

frequency  increment  of  wave  component  i 

phase  lead  of  motion  RAO  for  mode  j 

phase  lead  of  incident  wave  in  translating  earth  coordinate  system 

phase  lead  of  incident  wave  in  fixed  coordinate  system 

phase  lead  of  incident  wave  component  i 

incident  z  location  of  wave  surface 

motion  displacement  in  mode  j 

complex  motion  RAO  for  mode  j 

wave  heading  (to)  in  fixed  coordinate  system 

local  displacement  in  direction  j 

standard  deviation 

ship  heading  (toward)  in  fixed  coordinate  system 

incident  wave  frequency 

incident  frequency  of  wave  component  i 
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wave  encounter  frequency 
peak  wave  frequency 
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Annex  A 

ASCII  Post-Processing  Output  File  Format 
from  SHIPM07 


The  most  recent  version  (7.3)  of  DREA’s  SHIPMO  program  [1]  writes  output  post¬ 
processor  files  in  ASCII  rather  than  binary  format.  ASCII  files  are  compatible 
across  platforms  and  can  be  read  by  programs  written  in  any  computer  language. 
SHIPMO  version  7.3  uses  Fortran  90  allocatable  arrays,  which  greatly  reduces 
the  size  of  output  post-processor  files.  For  SHIPMO  7.3  to  write  an  ASCII  post¬ 
processor  file,  control  variable  OUTPPR  of  Record  (b)  of  the  SHIPMO  input  file 
must  be  set  to  OUTASCII. 

Record  1  -  Source  Program  Name 

Format  (A20) 

PROGRAM  Ch*20  Character  string  “SHIPM07.3”  identifying  program  name. 


Record  2  -  Title 

Format  (A  100) 

TITFE  Ch*100  Run  title  with  time  and  date  appended. 


Record  3  -  Control  Variables 

Format  (A20,  IX,  A20,  IX,  A20,  IX,  A20) 


OUTSYS 

Ch*20 

WATERTYPE 

Ch*20 

HYMETHOD 

Ch*20 

SPECTRUM 

Ch*20 

Output  units  (BRITISH  or  METRIC). 


Water  type  (SAFT  or  FRESH). 


Method  for  computing  hydrodynamic  coefficients  (BOUND- 
2D  or  CONFORM AF). 

Spectrum  type  (REGUFAR,  QUADRATIC,  BRETSCHNEI- 
DER,  JONSWAP,  OCHIHUBBFE,  INPUTUNI,  TENPARA- 
METER,  ODGPHINDCAST,  or  INPUTDIR) 
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Record  4  -  Ship  Dimensions 

Format  (F12.3,  IX,  14,  6(1X,  F12.3)) 


SHIPLEN 

Real 

Ship  length  between  perpendiculars  (ft  or  m). 

NST 

Integer 

Number  of  stations  (typically  21). 

DISP 

Real 

Ship  displacement  (tons  or  tonnes). 

LCBFP 

Real 

Distance  from  forward  perpendicular  to  LCB  (ft  or  m). 

DRAFTMID 

Real 

Draft  at  midships  (ft  or  m). 

TRIMS  T 

Real 

Trim  by  stem  (ft  or  m). 

KG 

Real 

Height  of  CG  above  baseline  (ft  or  m). 

BEAMMID 

Real 

Beam  at  midships  (ft  or  m). 

Record  5  -  Physical  Constants 

Format  (F12.3,  IX,  F12.3) 

WDENSITY  Real  Water  density  (slugs/ft3  or  kg/m3). 

GRAVITY  Real  Gravitational  acceleration  (ft/s2  or  m/s2). 

Record  6  -  Number  of  Stations  for  Sea  Load  Computations 

Format  (14) 

NSTLOAD  Integer  Number  of  stations  for  sea  load  computations. 

Record  7  -  Stations  for  Sea  Load  Computations 

Written  only  if  NSTLOAD  >  1  in  Record  (6) 

Format  (10(F10.3,  IX)) 

XSTLOAD(I)  Real  Array  of  dimension  NSTLOAD  with  station  numbers  for 

sea  loads. 
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Record  8  -  Number  of  Positions  for  Seakeeping  Calculations 

Format  (14) 

NPOS  Integer  Number  of  positions  for  seakeeping  computations. 


Record  9  -  Position  Locations  for  Seakeeping  Calculations 

Written  only  if  NPOS  >  1  in  Record  (8) 

Format  (9(F10.3,  IX)) 


XSTPOS(I) 

Real 

Array  of  dimension  NPOS  with  station  numbers  of  seakeep¬ 
ing  positions. 

YPOS(I) 

Real 

Array  of  dimension  NPOS  with  horizontal  coordinates  of 
seakeeping  positions  (ft  or  m). 

ZPOSBL(I) 

Real 

Array  of  dimension  NPOS  with  vertical  coordinates  of  sea¬ 
keeping  positions  relative  to  baseline  (ft  or  m). 

Record  10  -  Number  of  Ship  Speeds 

Format  (14) 

NSPEED  Integer  Number  of  ship  speeds. 


Record  11  -  Ship  Speeds 

Format  (10(F10.3,  IX)) 

SPDKNOT(I)  Real  Array  of  dimension  NSPEED  with  ship  speeds  (knots). 


Record  12  -  Number  of  Sea  Directions 

Format  (14) 

NSEADIR  Integer  Number  of  sea  directions. 

Record  13  -  Sea  Directions 

Format  (10(F10.3,  IX)) 

SEADIR(I)  Real  Array  of  dimension  NSEADIR  with  sea  directions. 
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Record  14  -  Number  of  Seaways 

Format  (14) 

NSEAWAY  Integer  Number  of  seaways. 


Record  15  -  Seaways 

Format  (10(F10.3,  IX)) 

HSW(I)  Real  Array  of  dimension  NSEAWAY  with  significant  wave  heights 

(ft  or  m). 

TSW(I)  Real  Array  of  dimension  NSEAWAY  with  characteristic  wave  pe¬ 
riods  (s). 


Record  16  -  Wave  frequencies 

Format  (14,  IX,  F9.3,  IX,  F9.3) 

NWVFREQ  Integer  Number  of  wave  frequencies. 

WVFREQMIN  Real  Minimum  wave  frequency  (rad/s). 

DWVFREQ  Real  Wave  frequency  increment  (rad/s). 

Ship  Motions  and  Loads  in  Regular  and  Irregular  Seas  for  each 
Ship  Speed 

The  remaining  sequence  of  records  is  repeated  for  each  ship  speed.  For  regular 
seas  and  irregular  seas  with  a  principal  wave  direction,  the  sea  direction  index  J 
in  the  motion  and  load  arrays  corresponds  to  the  sea  direction  relative  to  the  ship 
forward  speed.  For  irregular  seas  without  a  principal  wave  direction  (SPECTRUM 
=  TENPARAMETER,  ODGPHINDCAST,  or  INPUTDIR),  the  direction  index  J 
corresponds  to  the  ship  compass  heading  for  the  irregular  responses  only. 

Record  17  -  Label  for  Ship  Speed 

Format  (A40) 

SPEEDLABEL  Ch*40  “Results  for  Following  Ship  Speed  (knots)” 
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Record  18  -  Ship  Speed 

Format  (F9.3) 

SPEEDKT  Real  Ship  speed  (knots) 


Record  19  -  Label  for  Complex  Motion  RAOs  in  Regular  Waves 

Format  (A3  6) 

MOTIONRAOL AB EL  Ch*36  “Complex  Motion  RAOs  in  Regular  Waves” 


Record  20  -  Complex  Motions  in  Unit  Amplitude  Regular  Waves  for  Seaway  1 

Format  (5(’(\  E12.6,  E12.6,  ’)’)) 


CSURGEREG(U) 

Cmplx 

CSWAYREG(I,J) 

Cmplx 

CHEAVEREG(U) 

Cmplx 

CROLLREG(I,J) 

Cmplx 

CPITCHREG(U) 

Cmplx 

CYAWREG(U) 

Cmplx 

CRUDREG(I,J) 

Cmplx 

CSTABREG(IJ) 

Cmplx 

Complex  surge  amplitude  (ft/ft  or  m/m). 

Complex  sway  amplitude  (ft/ft  or  m/m). 

Complex  heave  amplitude  (ft/ft  or  m/m). 

Complex  roll  amplitude  (rad/ft  or  rad/m). 

Complex  pitch  amplitude  (rad/ft  or  rad/m). 

Complex  yaw  amplitude  (rad/ft  or  rad/m). 

Complex  rudder  amplitude  (rad/ft  or  rad/m). 

Complex  fin  or  U-tube  tank  amplitude  (rad/ft  or  rad/m). 


Notes:  The  above  arrays  have  dimensions  (NWVFREQ,  NSEADIR). 

Regular  wave  motions  and  loads  are  given  only  for  the  first  seaway 
regardless  of  the  number  of  seaways. 


Record  21  -  Label  for  Complex  Sea  Loads  in  Regular  Waves 
Written  only  if  NSTLOAD  >  1  in  Record  6 
Format  (A3 8) 

LOADRAOLABEL  Ch*38  “Complex  Sea  Load  RAOs  in  Regular  Waves” 
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Record  22  -  Complex  Sea  Loads  in  Unit  Amplitude  Regular  Waves  for  Seaway 

1 

Written  only  if  NSTLOAD  geq  1  in  Record  6 
Format  (5(’(\  E12.6,  E12.6,  IX)) 


CHORSHRREG(I,J,K) 

Cmplx 

CVERTSHRREG(I,J,K) 

Cmplx 

CTORSIONREG(I,J,K) 

Cmplx 

CVERTBNDREG(I,J,K) 

Cmplx 

CHORBNDREG(I,J,K) 

Cmplx 

Complex  horizontal  shear  amplitude  (lb/ft  or  N 

/m). 

Complex  vertical  shear  amplitude  (lb/ft  or  N 

/m). 

Complex  torsional  amplitude  (lb- ft/ft  or  N-m/m). 

Complex  vertical  bending  moment  amplitude 
(lb- ft/ft  or  N-m/m). 

Complex  horizontal  bending  moment  amplitude 
(lb- ft/ft  or  N-m/m). 


Note:  The  above  arrays  have  dimensions  (NWVFREQ,  NSEADIR,  NST¬ 

LOAD). 


Within  the  motions  and  loads  written  for  each  ship  speed,  the  following  sequence 
of  RMS  responses  in  irregular  seas  is  repeated  for  each  seaway  if  SPECTRUM  is 
not  equal  to  REGULAR. 

Record  23  -  Label  for  RMS  Motions  at  CG  in  Irregular  Waves 

Format  (A3  6) 

RMSCGLABEL  Ch*36  “RMS  Motions  at  CG  in  Irregular  Waves” 


Record  24  -  RMS  Motions  at  Ship  Centre  of  Gravity 

Format  (10(E12.6,  IX)) 


RMSSURGE(I) 

Real 

RMS  surge  (ft  or  m). 

RMSSWAY(I) 

Real 

RMS  sway  (ft  or  m). 

RMSHEAVE(I) 

Real 

RMS  heave  (ft  or  m). 

RMSROLL(I) 

Real 

RMS  roll  (degrees). 
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RMSPITCH(I) 

Real 

RMS  pitch  (degrees). 

RMSYAW(I) 

Real 

RMS  yaw  (degrees). 

RMSRUD(I) 

Real 

RMS  rudder  deflection  (degrees). 

RMS  STAB  (I) 

Real 

RMS  stabilizer  deflection  (degrees). 

RMS  S  URGE  VEL(I) 

Real 

RMS  surge  velocity  (ft/s  or  m/s). 

RMSSWAYVEL(I) 

Real 

RMS  sway  velocity  (ft/s  or  m/s). 

RMS  HE  AVE  VEL(I) 

Real 

RMS  heave  velocity  (ft/s  or  m/s). 

RMSROLLVEL(I) 

Real 

RMS  roll  velocity  (deg/s). 

RMSPITCHVEL(I) 

Real 

RMS  pitch  velocity  (deg/s). 

RMSYAWVEL(I) 

Real 

RMS  yaw  velocity  (deg/s). 

RMSRUDVEL(I) 

Real 

RMS  rudder  velocity  (deg/s). 

RMSSTABVEL(I) 

Real 

RMS  stabilizer  velocity  (deg/s). 

RMSSURGEACC(I) 

Real 

RMS  surge  acceleration  (g). 

RMS  S  WAYACC(I) 

Real 

RMS  sway  acceleration  (g). 

RMS  HE  AVE  ACC  (I) 

Real 

RMS  heave  acceleration  (g). 

RMS  ROLL  ACC  (I) 

Real 

RMS  roll  acceleration  (deg/s2). 

RMSPITCHACC(I) 

Real 

RMS  pitch  acceleration  (deg/s2). 

RMSYAWACC(I) 

Real 

RMS  yaw  acceleration  (deg/s2). 

RMSRUDACC(I) 

Real 

RMS  rudder  acceleration  (deg/s2). 

RMS  STAB  ACC  (I) 

Real 

RMS  stabilizer  acceleration  (deg/s2). 

Note:  The  above  arrays  have  dimension  (NSEADIR). 


Record  25  -  Label  for  RMS  Motions  at  Seakeeping  Positions  in  Irregular  Waves 

Written  only  if  NPOS  >  1  in  Record  (8) 

Format  (A54) 

RMSPOSLABEL  Ch*54  “RMS  Motions  at  Seakeeping  Positions  in  Irregular  Waves” 
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Record  26  -  RMS  Motions  at  Seakeeping  Positions 

Written  only  if  NPOS  >  1  in  Record  (8) 


Format  (10(E12.6,  IX)) 

RMSVERTDISP(I,J) 

Real 

RMS  vertical  displacement  (ft  or  m). 

RMSVERTVEL(U) 

Real 

RMS  vertical  velocity  (ft/s  or  m/s). 

RMSVERTACC(I,J) 

Real 

RMS  vertical  acceleration  (g). 

RMSRELDISP(I,J) 

Real 

RMS  relative  vertical  displacement  (ft  or  m). 

RMSRELVEL(I,J) 

Real 

RMS  relative  vertical  velocity  (ft/s  or  m/s). 

RMSLATDISP(I,J) 

Real 

RMS  lateral  displacement  (ft  or  m). 

RMSLATVEL(U) 

Real 

RMS  lateral  velocity  (ft/s  or  m/s). 

RMSLATACC(I,J) 

Real 

RMS  lateral  acceleration  (g). 

RMSLATFE(U) 

Real 

RMS  lateral  force  estimator  (g). 

RMSGFEPORT(IJ) 

Real 

RMS  force  estimator  function  for  tips  or  slides 
to  port  (g). 

RMSGFEPORTVEL(U) 

Real 

RMS  time  derivative  of  port  estimator  (g/s). 

RMS  GFES  TAR(I,  J) 

Real 

RMS  force  estimator  function  for  tips  or  slides 
to  starboard  (g). 

RMSGFESTARVEL(I,J) 

Real 

RMS  time  derivative  of  starboard  estimator  (g/s). 

RMSLONGDISP(I,J) 

Real 

RMS  longitudinal  displacement  (ft  or  m). 

RMSLONGVEL(I,J) 

Real 

RMS  longitudinal  velocity  (ft/s  or  m/s). 

RMSLONGACC(I,J) 

Real 

RMS  longitudinal  acceleration  (g). 

RMSLONGFE(IJ) 

Real 

RMS  longitudinal  force  estimator  (g). 

RMSGFEFORE(U) 

Real 

RMS  force  estimator  function  for  tips  or  slides 
forward  (g). 

RMS  GFEF  ORE  VEL(I,  J) 

Real 

RMS  time  derivative  of  forward  estimator  (g/s). 

RMS  GFEAFT  (I,  J) 

Real 

RMS  force  estimator  function  for  tips  or  slides 
aft  (g). 
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RMSGFEAFTVEL(I,J)  Real  RMS  time  derivative  of  aft  estimator  (g/s). 
Note:  The  above  arrays  have  dimensions  (NSEADIR,  NPOS). 


Record  27  -  Label  for  RMS  Sea  Loads  in  Irregular  Waves 

Written  only  if  NSTLOAD  >  1  in  Record  (6) 

Format  (A3  2) 


RMSLOADLABEL  Ch*32  “RMS  Sea  Loads  in  Irregular  Waves” 


Record  28  -  RMS  Sea  Loads 

Written  only  if  NSTLOAD  >  1  in  Record  (6) 
Format  (10(E12.6,  IX)) 


RMS  HORS  HR(I,J) 

Real 

RMS  horizontal  shear  (lb  or  N). 

RMSHORSHRVEL(IJ) 

Real 

RMS  time  derivative  of  horizontal  shear  (lb/s 
or  N/s). 

RMSVERTSHR(U) 

Real 

RMS  vertical  shear  (lb  or  N). 

RMS  VERT S  HRVEL(I,  J) 

Real 

RMS  time  derivative  of  vertical  shear  (lb/s  or 
N/s). 

RMSTORSION(IJ) 

Real 

RMS  torsion  (lb -ft  or  N-m). 

RMSTORSIONVEL(U) 

Real 

RMS  time  derivative  of  torsion  (lb- ft/s  or  N-m 
/s). 

RMSVERTBND(I,J) 

Real 

RMS  vertical  bending  moment  (lb -ft  or  N-m). 

RMS  VERTB  NDVEL(I,J) 

Real 

RMS  time  derivative  of  vertical  bending  mo¬ 
ment  (lb- ft/s  or  N-m/s). 

RMSHORBND(IJ) 

Real 

RMS  horizontal  bending  moment  (lb -ft  or  N-m). 

RMSHORBNDVEL(U) 

Real 

RMS  time  derivative  of  horizontal  bending  mo¬ 
ment  (lb- ft/s  or  N-m/s). 

Notes:  The  above  arrays  have  dimensions  (NSEADIR,  NSTLOAD). 
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Annex  B 

SHIPM07  Input  File  for  Frigate 


SHIPM07  input  file  of  frigate  for  ship  defence  models,  L  =  120  m 
METRIC  METRIC  SALT  NOSPEEDCOR  NOSWELLCOR  MOTION  OUTASCII  <-  Controls 


f rigateRao . dat 


NOLOAD  NORAW 
0.2  2.0  0.1 
NOSAVEHY  BOUND 2D 
0.2  6.0  0.2 
BRETSCHNEIDER 
13  0.0 

0  15  30  45  60  75 
1 

3.25  9.7 

7  0.0  5.0  10.0  15 
120.0  6.0  30.0 
GMCOMP  DRYROLLRG 
4.8 

OFFSETS 
21  21  1.0  1.0 
0.0 


LATLONG  HVCOR 


<- 


< - Additional  flags 

<  - Wave  frequencies 

< -  Hydrodynamic  options 

<  -  Encounter  frequencies 

<  - Wave  spectrum 

-  #  of  sea  dirs,  spread  angle 

90  105  120  135  150  165  180  < -  Sea  direction(s) 

< -  Number  of  seaways 

<  -  Wave  height,  period 

.0  20.0  25.0  30.0  < -  #  of  ship  speeds,  speeds 

<-  Ship  length,  KG,  pitch  gyradius 
< —  GM  and  roll  gyradius  flags 
< —  Dry  roll  radius  of  gyration 

<  - Hull  definition  flag 

<  -  #  of  stations,  scale  factors 
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0. 

.000 

0 

.030 

0. 

.036 

0, 

.048 

0. 

.065 

0, 

.090 

0. 

.122 

0. 

,161 

0. 

.208 

0. 

.262 

0 

.324 

0. 

.394 

0, 

.471 

0. 

.558 

0, 

.654 

0. 

.760 

0. 

,875 

0. 

.997 

1. 

.127 

1 

.264 

1. 

.409 

1, 

.560 

1. 

.718 

1. 

.881 

2. 

.049 

2. 

,300 

2. 

.512 

2. 

.725 

2 

.751 
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